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We investigate the role of a scalar tetraquark state for the description of nuclear matter within the 
parity doublet model in the mirror assignment. In the dilatation-invariant version of the model a 
nucleon-nucleon interaction term mediated by the lightest scalar tetraquark field naturally emerges. 
At nonzero density one has, beyond the usual chiral condensate, also a tetraquark condensate. The 
behavior of both condensates and the restoration of chiral symmetry at high density are studied. 
' It is shown that this additional scalar degree of freedom affects non negligibly the properties of the 

, medium. 
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I. INTRODUCTION 



The properties of strongly interacting matter at finite baryon density have been widely investigated in the past by 
use of chiral models [H, 0- This paper aims to study the effect of a light tetraquark field on such system. This subject 
is interesting for two reasons: (i) many works on light mesons spectroscopy show that the light scalar resonances below 
I , 1 GeV listed in the PDG [|[ can be successfully explained as a nonet of tetraquark states, see Refs. 0-0] and refs. 
Q-r therein, (ii) The nucleon can be modeled as a bound state of a quark and a (good) diquark S| . It is then conceivable 
that two nucleons, in addition to the usual quark exchange leading to a quark-antiquark meson as intermediate state, 
interact via an exchange of a diquark, which leads to a tetraquark as mediator, see Fig. [T] for a pictorial representation. 

In order to perform this study we use the linear sigma model for Nj = 2 described in Ref. @, in which the nucleon N 
and its chiral partner N* form a baryon doublet, (N, N*), where N* is usually identified with the resonance iV"(1535) 
0. The doublet is introduced in the so-called mirror assignment, first discussed in Ref. and extensively analyzed 
in Refs. foL ITTI. ll2T| . The particularity of the mirror assignment is the possibility to introduce a chirally invariant mass 
term ~ mo, which does not originate from the quark condensate. In the framework of dilatation invariant interactions, 
I/-) I this term originates from the condensation of two further scalar-isoscalar states: the dilaton/glueball field and the 
tetraquark field. In the present work we neglect the effect of the glueball since, due to its relatively high mass of 
about 1.5 GeV [l3|, its exchange between nucleons is negligible in a first approximation. On the contrary, a light 
tetraquark state with a mass of about 600 MeV and identified predominantly with the resonance /o(600) is potentially 
very interesting for the properties of nuclear matter. In addition, the usual (pseudo)scalar and (axial-)vector quark- 
antiquark mesons are present in the model. 

Applications of the parity doublet model to nuclear matter and neutron stars were studied in Refs. [rj - [l7j were 
it was shown that, at variance with the normal linear sigma model, it is able to describe the saturation of nuclear 
matter and it predicts a maximum mass for neutron stars compatible with observations. Some problems however arise 
from these studies: the nuclear matter compressibility turns out to be larger than the measured one and the value 
of mo needed to describe saturation is large, ~ 800 MeV, and thus in disagreement with the results of the analysis 
of Ref. [9j in which a fit of vacuum properties resulted in mo ~ 500 MeV. Moreover the mass of the scalar-isoscalar 
meson responsible for the nucleon-nucleon attraction turns out to be very small, m a ~ 350 MeV, a value that does 
not correspond to any particle in the PDG. 

In the present study, we remove these inconsistencies of the parity doublet model between vacuum physics and finite 
density physics by regarding the resonance /o(600) as a predominantly tetraquark state and the resonance /o(1370) as 
the chiral partner of the pion, hence a predominantly quark-antiquark state. We show that nuclear matter saturation 
can be correctly reproduced with values for the compressibility compatible with data. The value of mo needed for 
fitting high density properties is the same of the one indicated by Ref. @ in vacuum decays studies i.e. mo ~ 500 
MeV. Quite remarkably, the need of having two scalars within the parity doublet model is in agreement with the 
implementation, in nucleon-nucleon potentials as the one of Bonn |18| , of two scalar mesons having masses of ~ 500 
and - 1200 MeV. 

The paper is organized as follows: in Sec. II we introduce the Lagrangian of the model, in Sec. Ill we obtain 
the corresponding thermodynamic potential in the mean field approximation, in Sec. IV we presents our results for 
nuclear matter and, finally, in Sec. V we draw our conclusions. 
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FIG. 1: The left picture represents the exchange of quarks between nucleons leading the nucleon-nucleon interactions mediated 
by quark-antiquark mesons. The right picture represents the exchange of diquarks leading to an additional contribution to the 
nucleon-nucleon potential mediated by tetraquark mesonic states. 



Our units are h = c = 1, the metric tensor is g^ v = diag(+, — , — , — ). 



II. THE PARITY DOUBLET MODEL IN THE BARYON SECTOR 



We present here the chirally symmetric linear sigma model with scalar, pseudoscalar, vector, axial- vector mesons, 
the nucleon and its chiral partner Q. The scalar and pseudoscalar fields are included in the matrix 

3 

$ = ^2 <f>ata = (<r + ir/N) t + (So + in) ■ t, (1) 

a=0 

where t = r/2, with the vector of Pauli matrices r, and t° = I2/2. Under the global U(2)r x U(2)l chiral symmetry, 
$ transforms as $ — > Ul$>U R) where Ul and Ur are 2x2 unitary matrices. The vector and axial-vector fields are 
represented by the matrices 

3 3 

v = Y d = + ? ■ A^ = j2^t a = /r t Q + oi" . r . (2) 

a=0 a=0 

From these fields we define right- and left-handed vector fields = — A 11 , L M = + A^. Under global 
U(2) R x U(2) L transformations, these fields transform as R» -> UrR^U^ , -> U L L»ul. 

The identification of mesons with particles listed in the PDG 3] is as follows: the fields n and t]n correspond to the 
pion and the SU(2) counterpart of the rj meson, t]n = {uu + dd)/y/2, with a mass of about 700 MeV. This value can 
be obtained by "unmixing" the physical r\ and rf mesons, which also contain ss contributions. The vector fields co 11 
and f?* represent the resonances w(782) and p(770) and the axial-vector fields f£ and di^ represent the resonances 
/i(1285) and ai(1260). Two possibilities for the identification of the a and do fields exist: {/o(600), ao(980)} and 
{/o(1370), a (1450)}. The first assignment is however unfavoured 0, H3] (for a general discussion of the issue of 
scalar mesons see also Refs. [2l[ and refs. therein), while the second is in agreement with the phenomenology. 

The Lagrangian describing the meson fields is presented in detail in Ref. 0, For our purposes we notice that: 

(i) The chiral condensate <7o = (0 \a\ 0) = Z f n emerges upon spontaneous chiral symmetry breaking in the mesonic 
sector. The parameter /„• = 92.4 MeV is the pion decay constant and Z is the wave-function renormalization constant 
of the pseudoscalar fields [l9|, [22| and takes the value Z = 1.67 ± 0.2, which can be fixed from the process a\ — > nj. 

(ii) In the case Nf = 2 only one tetraquark state x exists. It can be coupled to the model following Refs. @, H^. (iii) 
The dilaton/glucball field G can also be easily added to the meson sector according to the requirement of dilatation 
invariance and the corresponding trace anomaly [2(|. (For a general discussion of the glueball see Refs. [2l|, [24| and 
refs. therein.) 

We now turn to the baryon sector, in which we have the baryon doublets "Fi and ^2, where \&i has positive parity 
and negative parity. In the mirror assignment they transform as follows: 
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FIG. 2: Energy per baryon as a function of the density for the case mo = 500 MeV for the chiral symmetry broken solution 
(solid line) and the chiral symmetry restored solution (dashed line). The dotted line represent the Maxwell construction which 
matches the two solutions. The model can describe the saturation of nuclear matter. 



> UrVir , *il > Ul^IL , *2fl > U L t>2R , *2L > U R ^ 2 L , (3) 

i.e., ^2 transforms in a "mirror way" under chiral transformations (TollTlj. These field transformations allow to write 
down a chiral baryonic Lagrangian: 

£bar = *lii7/^ll*l£ + ^lRH^DlR^lR + *2L*7pi^B*2L + ^ 2R.il 2 R 

- 91 + ^IR&^IL) - ?2 (*2L$ t *2_R + *2fl$*2L) 

- (flX + bG)(^ 1L ^2R - *17?*2L - *2L*li? + *2_R*li) , (4) 

where D^ R = <9 M — iciR^, D^ L = <9 M — iciL^, and Dt^ R = — ic 2 Rf ', D^ L = — ic 2 ^ are the covariant derivatives 
for the nucleonic fields, with the coupling constants c\ and C2- The interactions of the baryonic fields with the scalar 
and pseudoscalar mesons are parametrized by <?i and 32 ■ The last term in Eq. ((4]) generates a mass term when the 
tetraquark field % an d the glueball field G condense: 

m = axvac + bG vac ■ (5) 

In Ref. [9| the quantity too has been obtained through a fit procedure to known experimental and lattice quantities, 
obtaining: 

m = 460 ± 130 MeV . (6) 

In this paper we work under the simplified assumption 6 = 0, i.e. the parameter too is saturated by the tetraquark 
condensate and the corresponding nucleon-nucleon interaction is depicted in Fig. [TJ right side. In this context the 
coupling constant a is fixed as soon as the tetraquark condensate is specified, see next section. The physical fields N 
and N* are related to the spinors ^1 and ^2 through: 

*i = J-_ {Ne 5 ' 2 + l5 N*e^ 2 ) , * 2 = J__ ( l5 Ne- 8 ' 2 - N*e s ^ 2 ) . (7) 
V2coshd v2coshd 

The masses of the nucleons are obtained by diagonalizing the corresponding mass matrix in the Lagrangian. As 
function of the two condensates <r vac and Xvac they read: 




m N<N * = \l ( ' I ul ac + (aXvac) 2 ± 91 - ^ V vac- (8) 
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FIG. 3: Mean a (solid line) and \ (dashed line) fields. A strong first order phase transition is present at fi ~ 1 GeV. The 
transition point is the same for the two fields. Note, the minimal value for the chemical potential corresponds to nuclear matter. 



III. MEAN FIELD APPROXIMATION 

We now turn to the effective model in the mean field approximation suitable for describing nuclear matter. We 
thus retain only those fields which are relevant for such a study: the scalar fields a and x an d the vector field w. The 
effective Lagrangian is given by: 



1 ,22 , 1 2, 



1 



2 , 2 



+ -ma + -m^oT - -m x x 



A 



2 



-a +gxcr 



ea 



- ff^iiiT^*! - ^'^V^a + ax(^275^i - ^175*2) • (9) 
If we take their vacuum expectation values, a — > u + a vac and x ~^ X + Xvac, the potential for the fields a and x reads: 



-250-1,0 



A non diagonal term, which is proportional to the parameter g and mixes the quarkonium field a and the tetraquark 
field x, has emerged [30| . The physical fields, denoted as /i and s, are obtained through a standard diagonalization: 



h\ _ ( cos6» sin6» \ fx 
s J I — sin^o cos^o / V cr 

Then the physical masses of the scalar states are: 



1 4gcr„ a , 
in = — arctan 
2 



m 2 , — m 2 . 



to 2 = m 2 cos 2 6> + m 2 sin 2 6>o - 2ga vac sm(29 ) , 
to 2 = m 2 cos 2 + m?.sin flo + 2ga vac sm(29 ). 



(10) 



We identify here the predominately tetraquark state h with / (600) (the vacuum mass = 600 MeV is used) and 
the predominately quarkonium state s with / (1370) (the vacuum mass m s — 1300 MeV is used) . We now turn to 
the thermodynamic potential in the mean field approximation: 
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FIG. 4: Masses of the nucleons (solid line) and their chiral partners (dashed line) as function of the chemical potential. 



where the mesons term reads 
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-ma er 
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+ 9X& 



2 2 

2 m x x ■ 



The index i — N, N* denotes the nucleon type (positive and negative parity nucleons) , ji — 2 x 2 is the fermionic 
degeneracy (spin and isopsin), pp t are the Fermi momenta, E*(p) = \Jy\ + m^ 2 , and /i* = /x, — .g^wo : 



The single particle energy of each parity partner i is given by Ei{p) = E*(p) + g u L0Q. The mean mesons fields are 
obtained by minimizing the thermodynamic potential, i.e. by solving the following system of non-linear equations: 



d(Q/V) 
da 



dm 



\a 6 - m 2 a - e~ 2gxa + p* N (a, w , x) r ^ + P*n> w o , X) 



dm 



N 



da 



- m^uo - 9uj,nPn{o, uj , X) ~ 9^,n-Pn* (cr, u , X) = 



duj 
d(n/V) 



9 2 * i \ dm N dm* N 

ga +m x + PN\, a ^o,X)^ — + Pn'^^o^X)— — =0 • 



The scalar densities p* and the baryon densities pi are given by the expressions: 



7i 



d 3 p 



m 2 7T 2 



PFiEp. — mf In 



(11) 



(12) 



d 3 p 
(2^)3 



The parameters of the model are given by [lSj [23[ : 



2(^)2 



Z 2 



2g 2 



Z 2 



Z ' 



mom 2 . 



(13) 



(14) 



Numerically one has m^ = 139 MeV and /„• = 92.4 MeV. The parameter Z is the wave function renormalization 
of the pseudoscalar fields 0, l22 |. It is the only 'remnant' of the (axial-) vector mesons of the model. However, 
it should be stressed that the value of Z = 1.67 strongly affects the properties of nuclear matter, for instance 
reducing the value of the chiral symmetry explicit breaking term, as we will show in the following. As already 
discussed in Refs. ;lj|, the presence of (axial-)vector degrees of freedom in a chiral framework has non-negligible 
consequences on the (pseudo)scalar sector as well. The coupling constants of the nucleons and the scalar a meson are 
uniquely determined by the vacuum properties, m = 460 ± 130 MeV, mjy vac = 938 MeV and mN* ac = 1535 MeV: 
gi = 11.0 ± 1.5, <?2 = 18.8 ± 2.4 9] . Finally, g parametrizes the interactions of the tetraquark field with pions and the 
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FIG. 5: Compressibility as a function of mo (solid line) and allowed range for K as obtained by phenomenology (dashed lines). 
The constraint on the compressibility, together with vacuum decays properties, fixes mo to be in the range 500-600 MeV. The 
result (460 ± 130) MeV is shown (gray (yellow online) band). 



parameters g£ , 9u describe the interactions of the lo meson with the nucleons. For sake of simplicity we assume 



that g^ = #£ 2) = g u and thus, g^^ = g^,N* 
nuclear matter at saturation: 



The quantities g and g w are obtained from the following constraints of 



and 



=PF 



E 
A 



m N vaa 



-16 MeV. 



(15) 



where E/A is the energy per nucleon and pp N is the Fermi momentum of N. At the saturation point the latter equals 
value pf n — pf — 258 MeV. 



IV. RESULTS 



Let us fix mo to the intermediate value mo = 500 MeV. From Eqs. (10) and (15) we can determine the four free 
parameters of the model which turn out to be: g^ = 4.87, g — 450 MeV, m a — 1294 MeV, m x — 612 MeV (the latter 
two quantities generate the previously mentioned massed rrih = 600 MeV and m s = 1300 MeV). Quite remarkably 
the values of the masses of the two scalar mesons of our model are very close the values obtained within the "Bonn 
parametrization" of the nucleon-nucleon potential extracted by scattering data 18]. In Fig. [21 we show the energy 
per baryon E/A as a function of the baryon density p: one can notice the typical feature of the saturation of nuclear 
matter as a minimum of E/A located at p = po = 0.16 fm -3 and with a nucleon binding energy of 16 MeV. 

The value of the compressibility, computed through the derivative of the pressure P with respect to the density: 

= 194 MeV , (16) 

P=P0 

a value which is very close to the standard range indicated by phenomenology K — 200-300 MeV [25|, [26| . In the 
figure we also display the energy per baryon for the chiral symmetry restored solution of the mean field equations 
(dashed line) and the interpolating Maxwell construction (dotted line) which we will explain in the following. 

At high densities, our model predicts a strong first order chiral phase transition occurring at a chemical potential 
fi ~ 1 GeV as displayed in Fig. [3] where we show the mean scalar mesons fields (cr with solid line and \ with dashed 
line). The big jump in both the a and \ fields is computed by using a Maxwell construction: at a fixed value of the 
chemical potential the mean mesons fields are the solutions of Eqs. (11) which minimize the thermodynamic potential. 
In the energy per baryon - density diagram, the Maxwell construction corresponds to an hyperbolic branch (see Fig. 2, 
dotted line) given by the thermodynamic equation relating the pressure and the energy density: e/p = p cr it —Pcrit/p 
where p cr u and p cr it are the critical chemical potential and pressure. Notice that in absence of the tetraquark field 
and without the corrections related to Z [TEl ]. the phase transition, while being still a first order, is much weaker. 
The parameter Z = 1.67 indeed suppresses the chiral symmetry explicit breaking term and enhances the vacuum 
expectation value of a, thus making the results for the chiral phase transition similar to the ones obtained in the 
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FIG. 6: Critical densities associated with the chiral phase transition as functions of mo. pi (solid line) corresponds to the onset 
of the phase transition and p2 (dashed line) to the end (the densities are referred to the saturation density po). The first order 
phase transition becomes weaker for large values of mo. 

chiral limit. In our model the chiral phase transition is accompanied by the appearance of the chiral partners of 
the nucleon. While at low density only the nucleons are present, after the phase transitions both, the nucleons and 
their chiral partners, contribute the same amount. As one can notice in Fig. [4] in which we show the masses of the 
nucleons and their chiral partners as a function of the chemical potential, the nucleons become massless after the 
phase transition. This result is different from the one of Ref. [l5j in which the nucleon mass, at high density, saturates 
to the fixed value of tuq. On the other hand, in our model, mo itself is a function of the density and at the chiral 
phase transtion it vanishes since the tetraquark condensate vanishes (see Eq. 5). 

As a further step in our study we investigate the effect of varying the parameter mo. For each value of mo we 
compute the values of the four parameters described before: while we can describe saturation for the whole range of 
Eq. ©, the value of K turns out to be very small for small values of mo- We show the corresponding result in Fig. O 
together with the range of K indicated by phenomenology (dashed lines). By combing the information of Eq. (6) 
and the constraint obtained from the value of K we can further restrict the range of the allowed values for the "bare 
mass" of the nucleon 500 < m < 600 MeV. 

Finally, in Fig. [5J we show the values of the critical densities pi and p2 which correspond to the onset and the end of 
the chiral phase transition. One can notice that the strength of the phase transition decreases for large values of mo; 
within the allowed range of mo the density jump at the phase transition is anyway large : p%j p\ — 2-4. Remarkably, 
recent lattice results [22} show for a temperature T ~ 0.8T C a rather broad range of density for the mixed phase, 
Pil P\ ~ 2. Such studies are presently limited for temperatures T > 0.8T C , therefore a quantitative comparison with 
our results is premature. Nevertheless, the lattice simulations seem to indicate that the interval for the mixed increase 
further when diminishing the temperature T. Future detailed analysis on this issue is surely interesting. 

V. CONCLUSIONS 

In this work we have used the chiral model developed in Ref. [9] to study nuclear matter properties. The model 
contains, in addition to the standard quark-antiquark fields with (pseudo)scalar and (axial-) vector quantum numbers, 
also a light scalar-isoscalar tetraquark field. In the framework of the mirror assignment the latter can be easily 
coupled to the nucleon by using the general requirement of dilatation invariant interactions and contributes to the 
nucleon- nucleon interaction, see Fig . 1. 

Similarly to the results of Ref. [23] at finite temperature and vanishing density, we find that also at nonzero density 
a light tetraquark field has a strong influence on the medium properties of the system due to the interplay of two 
condensates, the tetraquark and the chiral (quark-antiquark) condensates. Interestingly, the described scenario and 
also the values of the coupling constants are in agreement with the work of Ref. [l8|, where two scalar fields are 
needed to describe nucleon- nucleon scattering data. As described in Ref. [28[, in such a scenario nuclear matter is a 
prerogative of our world with three colors (N c = 3), but would cease to exist as soon as a larger number of colors is 
considered. 

An important parameter of the model is mo, which describes the contribution to the nucleon mass which does not 
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stem from the chiral condensate and which, in the present study, is saturated by the tetraquark condensate. Its value 
has been fixed in the vacuum as m = 460 ±130 MeV. It is then remarkable that the compressibility is in agreement 
with the experiment in a compatible range of mo (see Fig. 5), thus showing that vacuum results and nonzero density 
properties can be understood within the same theoretical model. 

Further studies along the direction of the present work can be performed: (i) the contribution of the gluon condensate 
to too and of the dilaton/glueball field to the nucleon-nucleon interaction should be included. Although the glueball 
is heavier (~ 1.5 GeV) and should not affect the interaction of two nucleons, its inclusion is important in a theoretical 
framework which is based on the dilatation invariance of QCD and on its anomalous breaking, (ii) The study of 
asymmetric nuclear matter and its application to neutron stars is also important since, recently, very massive neutron 
stars have been discovered which can give useful constraints for the stiffness of the equation of state. Moreover, also 
the symmetry energy of nuclear matter at high density is a crucial quantity that we can investigate by introducing 
in our model the scalar isovector ao(980), as done in Ref. [29j], which in our assignment is also a tetraquark. (iii) The 
extension of the model to Nf = 3 is currently under investigation. The effects of nonzero density and temperature of 
the complete model with all the relevant degrees of freedoms represents an interesting outlook of the present work. 

The work of G.P. is supported by the Deutsche Forschungsgemeinschaft (DFG) under Grant No. PA 1780/2-1. Wc 
thank D.H. Rischke, J. Schaffner-Bielich and L. Bonanno for valuable discussions. 



[1] N.K. Glendenning, Compact Stars, 2nd ed. 2000, Springer- Verlag, New York, Berlin, Heidelberg. 

[2] I. Mishustin, J. Bondorf and M. Rho, Nucl. Phys. A 555 (1993) 215. P. Papazoglou, S. Schramm, J. Schaffner-Bielich, 
H. Stoecker, W. Greiner, Phys. Rev. C57 (1998) 2576-2588. P. Papazoglou, D. Zschiesche, S. Schramm, J. Schaffner- 
Bielich, H. Stoecker, W. Greiner, Phys. Rev. C59 (1999) 411-427. L. Bonanno, A. Drago, A. Lavagno, Phys. Rev. Lett. 
99 (2007) 242301. L. Bonanno, A. Drago, Phys. Rev. C79 (2009) 045801. E. K. Heide, S. Rudaz, P. J. Ellis, Nucl. 
Phys. A571 (1994) 713-732. 

[3] K. Nakamura et al. (Particle Data Group), J. Phys. G 37, 075021 (2010). 

[4] R. L. Jaffe, Phys. Rev. D15, 267 (1977). R. L. Jaffe, Phys. Rev. D15, 281 (1977). 

[5] L. Maiani, F. Piccinini, A. D. Polosa, and V. Riquer, Phys. Rev. Lett. 93, 212002 (2004), hep-ph/0407017. G. 't Hooft, 
G. Isidori, L. Maiani, A. Polosa, and V. Riqu er, Phys.Lett. B662, 424 (2008), 0801.2288. A. H. Fariborz, R. Jora and 
J. Schechter, Phys. Rev D 72 (2005) 034001 |arXiv:hep-ph/050 6170 . A. H. Fariborz, Int. J. Mod. Phys. A 19 (2004) 2095. 
[arXiv:hep-ph/0302133| . M. Napsuciale and S. Rodriguez, Phys. Rev. D 70 (2004) 094043. 

[6] F. Giacosa, Phys. Rev. D74, 014028 (2006), hep-ph/0605191. 

[7] F. Giacosa, Phys.Rev. D75, 054007 (2007), hep-ph/0611388. 

[8] G. Hellstern, R. Alkofer and H. Reinhardt, Nucl. Phys. A 625 ( 1997) 697 [arXiv:hep-ph/9706551| . G. Hellstern, R. Alkofer, 
M. Oettel and H. Reinhardt, Nucl. Phys. A 627 (1997) 679 |arXiv:hep'ph/9705267| . M. Oettel, R. Alkofer and L. von 
Smekal, Eur. Phys. J. A 8 (2000) 553 [arXiv:nucl-th/0006082| . A. Bender, C. D. Roberts and L. Von Smekal, Phys. 
Lett. B 38 (1996) 7 [arXiv:nucl-th/9602012| . P. Maris~and C. D. Roberts, Int. J. Mod. Phys. E 12 (2003) 297 
'arXiv:nucl-th/0301049| . U. Vogl and W. Weise, Prog. Part. Nucl. Phys. 27 (1991) 195. 



[9 
[10 

[11 

[12 

[13 

[14 
[15 
[16 

[17 
[18 
[19 



[20 
[21 

[22 
[23 
[24 



S. Gallas, F. Giacosa, and D. H. Rischke, Phys.Rev. D82, 014004 (2010), 0907.5084. 

B. W. Lee, Chiral Dynamics (Gordon and Breach, New York, 1972). 

C. E. Detar and T. Kunihiro, Phys. Rev. D39, 2805 (1989). 

D. Jido, M. Oka, and A. Hosaka, Prog.Theor.Phys. 106, 873 (2001), hep-ph/0110005. D. Jido, Y. Nemoto, M. Oka, and 
A. Hosaka, Nucl.Phys. A671, 471 (2000), hep-ph/9805306. 

Y. Chen, A. Alexandru, S. Dong, T. Draper, I. Horvath, et al., Phys.Rev. D73, 014516 (2006), hep-lat/05 10074. 
T. Hatsuda and M. Prakash, Phys. Lett. B224, 11 (1989). 

D. Zschiesche, L. Tolos, J. Schaffner-Bielich, and R. D. Pisarski, Phys.Rev. C75, 055202 (2007), nucl-th/0608044. 

V. Dexheimer, S. Schramm, and D. Zschiesche, Phys.Rev. C77, 025803 (2008), 0710.4192. V. Dexheimer, G. Pagliara, 

L. Tolos, J. Schaffner-Bielich, and S. Schramm, Eur.Phys.J. A38, 105 (2008), 0805.3301. 

C. Sasaki and I. Mishustin, Phys. Rev. C 82 (2010) 035204, 1005.4811. 
R. Machleidt, Phys.Rev. C63, 024001 (2001), nucl-th/0006014. 

D. Parganlija, F. Giacosa, and D. H. Rischke, Phys.Rev. D82, 054024 (2010), 1003.4934; D. Parganlija, F. Giacosa, 
and D. H. Rischke, AIP Conf.Proc. 1030, 160 (2008), 0804.3949; D. Parganlija, F. Giacosa, and D. H. Rischke, PoS 
CONFINEMENTS, 070 (2008), 0812.2183 [U 

S. Janowski, D. Parganlija, F. Giacosa and D. H. Rischke, arXiv:1103.3238 [hep-ph]. 

C. Amsler and N. A. Tornqvist, Phys. Rept. 389, 61 (2004). E. Kle mpt and A. Zait sev, Phys. Rept. 454 (2007) 1 
|arXiv:0708.4016l [hep-ph]]. F. Giacosa, Phys. Rev. D 80 (2009) 074028 |arXiv:0903.448T1 [hep-ph]'] . 
S. Struber and D. H. Rischke, Phys.Rev. D77, 085004 (2008), 0708.2389. 

A. Heinz, S. Struber, F. Giacosa, and D. H. Rischke, Phys.Rev. D79, 037502 (2009), 0805.1134. 

C. Amsler and F. E. Close, Phy s. Rev. D 53 (1996) 2 95 [arXiv:hep-ph /9507326| . W. J. Lee and D. Weingarten, 
Phys. Rev. D 61, 014015 (2000). [arXiv:hep-lat /9910008| ; F. E. Close and A. Kirk, Eur. Phys. J. C 21, 531 (2001). 



9 



arXiv:hep-ph/0103173 . F. Giacosa, T. Gutsche, V. E. Lyubovitskij and A. Faessler, Phys. Rev. D 72, 094006 (2005). 
arXiv:hep- ph/0509247 . F. Giacosa, T. Gutsche, V. E. Lyubovitskij and A. Faessler, Phys. L ett. B 622, 277 (2005 ) 
arXiv:hep-ph/0504033 . F. Giacosa, T. Gutsche and A. Faessler, Phys. Rev . C 71, 0252 02 (2005) [arXiv:hep-ph/040 8085 . 

.1. Y. Cheng, C. K. Chua and K. F. Liu, Phys. R ev. D 74 (2006) 094005 |arXiv:hep-ph/0607206| . Y. Chen et al, Phys. 

R ev. D 73 (2006) 014516. [arXiv:hep-lat /0510074] . V. Mathieu, N. Kochelev and V. Vento, Int. J. Mod. Phys. E 18 (2009) 

1 |arXiv:0810.4453l [hep-ph]]~ 
[25] D. H. Youngblood et al., Phys. Rev. C69, 034315 (2004). 

[26] C. Hartnack, J. Aichelin, and H. Oeschler, Int. J. Mod. Phys. D16, 153 (2007). 

[27] A. Li, A. Alexa ndru, K. F. Liu an d X. Meng, Phys. Rev. D 82 (2010) 054502 |arXiv: 1005.41581 [hep-lat]] . A. Li, A. Alexandra 
and K. F. Liu. larXiv: 1103.30451 [hep-ph]. A. Li, A. Alexandru, X. Meng and K. F. Liu [chi QCD Collaboration], Nucl. 
Phys. A 830 (2009) 633C arXiv:0908.1155 [hep-lat]]. 

[28] L. Bonanno and F. Giacosa, arXiv:1102.3367 [hep-ph]. 

[29] B. Liu, V. Greco, V. Baran, M. Colonna and M. Di Toro, Phys. Rev. C 65 (2002) 045201 [arXiv:nucl-th/0112034] . 

[30] The term responsible for the tetraquark-quarkonium interaction, and therefore also for the mixing, is gx& in Eq. @. It 
was introduced in Ref. Q because it is dominant in the large-iV c expansion. (Terms of the kind X 2(j2 would also possible, 
but they are here neglected because large- N c suppressed.) Note also that the term gx& 2 is the only non-dilatation invariant 
term in the chiral limit. In a fully dilatation invariant context one should write down the term aGx° 2 where G is the 
dilaton field. Upon condensation of the latter one obtains g = aGo- As long as the fluctuations of the dilaton are neglected 
it is justified to work with gx& 2 ■ 



